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Abstract 

We prove that it is NP-hard to decide whether two points in a polygonal domain with holes can 
be connected by a wire. This implies that finding any approximation to the shortest path for a long 
snake amidst polygonal obstacles is NP-hard. On the positive side, we show that snake's problem is 
"length-tractable": if the snake is "fat", i.e., its length/width ratio is small, the shortest path can be 
computed in polynomial time. 
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1 Introduction 

The most basic problem in VLSI and printed circuit 
board design is to connect two given points, s and t, 
by a shortest " thick" path avoiding a set of polygonal 
obstacles in the plane. The quarter-of-a-century-old 
approach to the problem is to inflate the obstacles 
by half the path width, and search for the shortest 
s-t path amidst the inflated obstacles [9]. The found 
path, when inflated, is the shortest thick s-t path. 

It went almost unnoticed that the thick path built 
by the above procedure may self-overlap (Fig. IT]): 
apart from our recent work on thick paths ^ , we only 
found one mention of the possibility of the overlap 
— Fig. 4 in [13]. (In a different context, Bereg and 
Kirkpatrick [5^, Fig. 2] also noted that Minkowski sum 
of a disk and a path may be not simply-connected.) 
When the path represents a thick wire connecting 
terminals on a VLSI chip or on a circuit board, self- 
overlap is undesirable as the wire must retain its width throughout. Thus, the objective in the basic wire 
routing problem should be to find the shortest non- self overlapping thick path. 

The problem shows up in other places as well. For instance, one may be interested in the optimal conveyor 
belt design: the belt is a non-selfoverlapping thick path. Our particular motivation comes from air traffic 
management where thick paths represent lanes for air traffic. Lane thickness equals to the minimum lateral 
separation standard, so that aircraft following different lanes stay sufficiently far apart to allow for errors in 
positioning and navigation. If an airlane self-overlaps, two aircraft following the lane may come too close to 
each other; thus it is desirable to find lanes without self-overlaps. 




Figure 1: Only a self-overlapping thick path exists. 
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1.1 Our contributions 

We prove a surprisingly strong negative result (Section [4]) : it is NP-hard even to decide whether there 
exists (possibly, arbitrarily long) s-t wire; this implies that no approximation to the shortest wire can be 
found in polynomial time (unless P=NP). 

Short Snakes Our intractability result means that in general it is NP-hard for a snake to wriggle its way 
amidst polygonal obstacles (assuming the snake is uncomfortable with squeezing itself). The good news for 
snakes is that in our hardness proof the sought wire is considerably long; i.e., the hardness of path finding 
applies only to long snakes. Our positive result (Section [3| is that for a bounded-length snake, the shortest 
path can be found in polynomial time (assuming real RAM and the ability to solve constant-size differential 
equations in constant time) by a Dijkstra-like traversal of the domain. 

1.2 Related work 

In VLSI numerous extensions and generalizations of the basic problem were considered. These include routing 
multiple paths, on several levels, and with different constraints and objectives. It is impossible to survey all 
literature on the subject; we will only mention the books [I9lf20] . 

In robotics thick paths were studied as routes for a circular robot. In this context, path self-overlap poses 
no problem as even a self-overlapping path may be traversed by the robot; that is, in contrast to VLSI, 
robotics research should not care about finding non-selfoverlapping paths. In [3], Chew gave an efficient 
algorithm for finding a shortest thick path in a polygonal domain. In a sense, our algorithm for shortest path 
for a short snake (Section^ may be viewed as an extension of Chew's. 

Motion planning for an object with few degrees of freedom may be approached with the cell decomposition 
techniques [TH [TB] . Closest to our bounded-length snake problem is the work on path planning for a segment 
(rod) (5, [61 [T71 [22] . Short snakes are also relevant to more recent applications of motor protein motion [TUl [23] . 

2 Snake Anatomy and Physiology 

In this section we introduce the notation and formulate our problem. 

Let P be an n-vertex polygonal domain with obstacles. For a planar set S let hdS denote the boundary 
of S, and for r > let < S >^^' denote the Minkowski sum of S with the radius-r open disk centered at 
the origin. Let P^ — P\ < hdP >^^' be P offset by r inside. The boundary of P*" consists of straight-line 
segments and arcs of circles of radius r centered on vertices of P. We call such (maximal) arcs r-slides. 

Let TT be a path within P^; let |7r| denote its length. A thick path H is the Minkowski sum 11 = < tt >(^\ 
The path tt is called the reference path of 11; the length of 11 is |7r|. 

A snake is a non-selfoverlapping thick path, i.e., a path which is a simply-connected region of the plane. 
The reference path of the snake is its spine (Fig. [2|. One of the endpoints of the spine is the snake's mouth 
m. The snake is a "rope" that "pulls itself by the head" : imagine that there are little legs (or a wheel, for a 
toy snake) located at to, by means of which the snake moves. The friction between the snake's body and the 
ground is high: any point p of the spine will move only when the path from the mouth to p is a "pulled-taut 
string", i.e., is a locally shortest path. That is, the snake always stays pulled taut against the obstacles (or 
itself). 

The input to our problem is the domain P, two points s, f ^ P^ ~ the "start" and the "food", a number 
L > - the length of the snake, and the initial direction of the snake at s. (Assume w.l.o.g. that s and / are 
at distance 1 from some vertex of P.). The goal is to find a path for the snake such that the snake's mouth 
starts at s and ends at /; the constraints are that the snake remains pulled taut and non-selfoverlapping 
throughout the motion. The objective is to minimize the distance traveled by the mouth, or equivalently, 
assuming constant speed of motion, the time until the snake reaches the food. 

Remark 1 . To simplify the formulation, we did not specify the initial configuration of the snake; a pedantic 
view could be to assume that the snake slithers in from a Riemann sheet glued along the diameter of < s >''^' 
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Figure 2: Spine, mouth, head, headcut, antimouth, tail, bends. A layer-2 tailed bend Bi does not have a 
point connected to v{Bi) with a length-2 segment fully lying within the snake. 

that is perpendicular to the initial direction of the snake's motion. Our results, both positive and negative, 
remain valid even if the initial configuration of the snake is part of the input. 

Remark 2. It is not true that all points of the spine necessarily follow the same path. 

Remark 3. Whoever guesses that the above model of a snake was developed just for FUN, is right. Nevertheless, 
the proposed problem formulation may be relevant also in more serious circumstances. It models, e.g., the 
path of a rope being pulled by its frontpoint through a polygonal domain. If it is a fire hose or a tube 
delivering life-saving medicine [3l [11], minimizing the time to reach a certain point seems like a natural 
objective (more important than, say, the work spent on pulling the tube). For another application, consider 
a chain of robots moving amidst obstacles. Each robot, except for the leader of the chain, has very simple 
program of following its predecessor - just keeping the distance to it. Then the robots form (approximately) 
a pulled taut thick string. 

Remark 4. It is possible to come up with problem instances in which no path for a pulled-taut snake exists, 
while a path for a snake that is not required to be taut, does. 

3 Shortest Path for a Fat Hippo 

When the snake is relaxed and its spine is a straight-line segment, the snake is the Minkowski sum of the 
length-L segment and the unit disk. Such sums are known as hippodromes [U [21 [71 [121 HI] j or hippos for short. 
We say that a hippo is fat if its length is constant: L — 0{1). In this section we show that for a fat hippo 
our problem can be solved in polynomial time. 



Overview of the approach 

Our algorithm is a Dijkstra-style search in an implicitly defined graph Q (Fig. [s]) : neither the nodes nor the 
edges of the graph are known in advance. Instead, Q is built incrementally, by propagating the labels from 




Figure 3: w is a node of G', pull the snake from v for distance L along "every possible" path. Endpoints of 
visibility edges (some of them are shown with solid circles) become nodes of G- 

the node v with the smallest temporary label (as in standard Dijkstra). The labels are propagated to nodes 
in the "L- visibility" region of v, which is what the snake "sees" while it slithers for distance L starting from 
V. In order to discover the nodes in the region, we pull the snake from v for length L along every possible 
combinatorial type of path; by a packing argument, there is only a small number of the types. The edges of G 
correspond to bitangents between the paths and 1-slides. We prove that the algorithm is polynomial-time by 
observing that the snake must travel for at least 27r — 2 before it "touches" itself with its head; this implies 
that the snake never "covers" any point of bdP with more than 0(L/7r) — 0(1) "layers", and hence there is 
a polynomial number of relevant bitangents. 

In the reminder of the section we elaborate on the algorithm's details. 

3.1 A Bit More Anatomy 

Consider the unit disk < m >'^^\ The part of the boundary of < m >(^) that is also the boundary of the 
snake is the unit semicircle whose diameter is perpendicular to the spine at m; we call the part the head and 
the diameter the headcut. The endpoint of the spine that is not the mouth is called the antimouth a. The 
part of the boundary of < a >^^^ that is also the boundary of the snake is called the tail. Refer to Fig. [2J 



3.2 Freeze! 

Let us have a closer look at the structure of the pulled-taut snake at any moment of time. Some pieces of 
the spine are straight-line segments. The segments are bitangents between the other, non-segment pieces 
supported by vertices of P, possibly via several "layers" of the snake. We call each such (maximal) piece B a 
bend; we denote the vertex that supports B by v{B), and say that v is responsible for B. 



Snake layers We show that each vertex v is responsible for 0(1) bends. Say that a bend B belongs to 
layer 1 if there exists a point b £ B such that |5t;(i?)| = 1. Recursively, B is layer-{k + 1) bend if it is not 
layer-fc and there exists a point b E B such that for some point b' at layer k, we have \b'b\ — 2 and b'b fully 
lies within the (closure of the) snake (Fig. ^. Let K be the maximum index of a layer. 

Lemma 3.1. K < [L/(27r - 2)] . 

Proof. Let bk,bk+i be points on layers k,k + 1 such that |6fcfefc+i| = 2. Let TT{bkbk+i) be the part of the 
spine between bk and bk+i- The Minkowski sum < TT{bkbk+i) > encloses at least one obstacle, h (or 
else the snake is not pulled taut). The perimeter of the inflated obstacle < h >^ ' is at least 2tt, thus 
\TT{bkbk+i)\ + \bk+ibk\ > 277. But \bk+ibk\ = 2, hence \Tr{bkbk+i)\ > 27r ~ 2. D 

As a corollary from the lemma, we have that life is very simple for a fat enough hippo: 

Corollary 3.2. For a > 27r — 2, an a-fat hippo can follow a shortest thick path without self-overlap. 

Tailed and headed bends One may wonder why we did not opt for a simpler definition of layer-fc bend B 
as one having a point b such that |6w(i?)| =2^ — 1 and bv{B) lies fully within the (closure of the) snake. The 
reason are two special kinds of bends which make the snake's portrait more complicated than in the case of 
an infinite-length snake. Specifically, we say that a bend B at layer k is tailed (resp. headed) ii < B >^ ~ ' 
touches the tail (resp. head). The simpler definition may not work for such bends (Fig. ^. 

Any bend that is not tailed or headed is an arc of a circle (actually, it is part of a slide). A tailed or 
headed bend i? is a different shape - string pulled taut against a ball touching v{B). 

Snake configuration The snake can be reconstructed in linear time as soon as the following is specified: 
(1) list of the vertices responsible for the bends; (2) for each bend, its layer; (3) for a headed bend - the vertex 
or the edge of P in contact with the head, and the slope of the headcut; (4) similar information for a tailed 
bend; (5) positions of the mouth m and antimouth a. We will call (1)~(5) the configuration of the snake. 

Lemma 3.3. The list (1) contains 0{n) vertices. 

Proof. Each vertex may be responsible for up to K bends, hence the total number of bends is 0{Kn) — 0{n). 
(In fact, since the part of the snake between consecutive bends leaning on one vertex, encloses at least one 
obstacle, every bend can be charged either to a vertex or to an obstacle; thus the number of bends is actually 
0{n + h) = 0{n), where h is the number of obstacles.) D 

3.3 Move! 

Suppose that we are given the configuration C of the snake at some time; let tt be the spine when the snake is 
in C and let mc be the position of the mouth in C. Suppose we are also given a path 7 for the mouth starting 
at niQ. Assume 7 has the following properties: (1) it is consistent with C in that the tangent to 7 at mc 
coincides with the tangent to n at mc', (2) it has a polynomial number of pieces, each of constant description 
complexity; (3) I7I < L. In what follows we assume that every path 7 has these properties. We claim that 
in polynomial time we can check whether 7 is a feasible path for the mouth, i.e., whether the snake stays 
obstacle- free when m is pulled along 7 (say, at unit speed). 

First we note that it is enough to check only whether the mouth moves feasibly. Indeed, the first time 
that the snake (possibly) becomes infeasible as m follows 7, the mouth is necessarily a part of "certificate of 
infeasibility" . This is so because the only way that the snake experiences "side pressure" is due to appearance 
of a headed bend. Other than that, any piece V of the snake is merely pulled by the preceding piece, V': 
either V exactly follows the same path that just was feasible for V' , or 7^ is a tailed bend B. Of course, the 
bend niorphs as the time passes, but only becomes "more feasible" with time, i.e., the free space around B 
increases - _B is pushed only by the tail, and the tail "moves away" with time. 

Let now 7' be a piece of 7. If we know how each piece V of the snake changes with time, we can test 
whether m stays feasible while following 7', by checking the feasibility against each V in turn. That is, 



while neither the configuration of the snake nor the piece 7' of 7 changes, the feasibility test can be done 
piece- versus-piece in constant time (assuming real RAM). Observe that overall there is only a polynomial 
number of configuration changes. Indeed, the configuration may change only due to one of the following events: 
(1) the tail starts to follow another feature of P, (2) a headed bend appears or changes its combinatorial 
structure, (3) a tailed bend disappears or changes its combinatorial structure. But each of the events (l)-(3) 
may happen only once per vertex-bend-layer triple; thus, by Lemmas |3.1| and |3.3| there is only a polynomial 
number of events. 

Tracking the configuration changes is easy given the way each piece changes with time. For event (1), we 
only have to know how the tail speed changes with the time between consecutive events (the tail speed is not 
necessarily constant, we elaborate on it in the next paragraphs). For event (2) we check what is the first time 
that the head collides with a piece or when a headed bend hits a vertex; all this can be done in polynomial 
time as there is only a linear number of candidate collisions. Event (3) is similar. The next event time is 
then the minimum of the event times over all the pieces. 

It remains to show how to determine the way each piece changes. Here the crucial role is played by 
the headed and tailed bends; again, it is the finiteness of the snake length that makes things involved. Let 
Bi — bib\ be the first such bend counting from m (Fig. [2]). Assume that Bi is a layer- A; tailed bend; the 
situation with a headed bend is actually simpler (because the mouth speed is constant). Let cibi and b'lc'i 
be the pieces adjacent to Bi - both are bitangents (straight-line segments, possibly of length) to Bi and 
adjacent bends. Every point of the spine between m and ci moves at speed 1, and none of the bends before 
Bi changes with time. 

To figure out what happens after ci, we have to solve a constant-size differential equation that describes 
the "propagation" of speed of motion of the spine. Specifically, let u(t) denote the speed at which the 
antimouth moves at time t. Knowing u{t) and knowing the initial position of the antimouth, we can write 
the antimouth position as a function of time, and hence we know < a >(2'=-2) g^g ^ function of time. The 
points 5i and b[ are points of tangency to < a >(^'^~^) from Ci and c'l; thus knowing < a >(2'=-2) .^^g know 
how the length |ci6i| -I- |7r(6i6'^)| -I- \b[c[\ changes with time. Knowing that, and recalling that at Ci the spine 
moves with unit speed, we can write what the spine speed at cf^ is as a function of r. 

Now, the spine speed does not change between c[ and the next point, C2, that is the start of the tangent 
to the next headed or tailed bend, -82- We perform at B2 the same operations as above, and get the speed of 
motion of the spine past the bend i?2- Continuing in this fashion, in the end, after going through all bends, 
we obtain some expression, £(u(t)) for the spine speed after the last bend. 

Finally, to close the loop, we solve the equation 

u{t)=£{u{t)) (1) 



Since there is only a constant number of layers (Lemma 3.1 1, there is only a constant number of the headed 
and tailed bends, and hence the expression f is a sum of a constant number of terms (each containing u{t) 
and Ju{T)dT), each of constant description complexity. In our computation model, we can solve the equation 
for u{t) in constant time. Substituting u{t) back into the formulae for the different bends, we obtain the 
spine as a function of time, as desired. 

3.4 See! 

Assume that at some point the snake is in configuration C. What happens next, as the snake follows the 
optimal path to /? Local optimality conditions imply that it will "wiggle around the obstacles" for some 
time and then "shoot" towards a vertex of P. We formalize this below. 

Final piece of anatomy The eye of the snake is collocated with the mouth. The snake can see a point 
X £ P^ if the segment mx lies fully within P^ and is tangent to the spine at m. Recall that P^ = P\< hdP >^^' 
where < bdP >^ Ms the Minkowski sum of bdP with open unit disk; hence mx can go along the boundary 
of P^ (with X being, e.g., an endpoint of a slide; see Fig. ^. The snake itself is transparent for its eye: we do 
not forbid mx to intersect the snake. 



Definition 3.4. A point p on bdP^ is L-visible from mc if there exists a path 7 for the mouth ending in a 
point m* such that m* sees p and m*p is tangent to bdP^ at p. We say that mcp is an L-visibility edge, or 
an L-edge for short. We say that 7 is the wiggle segment of mcP, and that 77i*p is the visibility segment of 



the edge. (We remind that we assume 7 enjoys the properties hsted in the beginning of Section 3.2 tangent 
at mc consistent with C, polynomial-size description, length < L.) 

To be on a (locally) optimal path, the mouth would like to follow an L-edge also past m* . This may not 
be feasible due to a conflict with the snake itself. However, such conflicts can be discovered "on-the-fly", 
as the mouth attempts to move along m*p. Specifically, try to move the mouth along m*p, as described in 
Section [33] If during the motion, the head collides with the snake, note what kind of bend the head collides 
with. If this bend B is not tailed, adjust the path for the mouth so that it is tangent to < -B >^^-', and follow 
the adjusted path. If by the time the mouth reaches < B >*-^% the bend B is "gone", i.e., the head "misses" 
the snake, we know that we are dealing with a tailed bend (or with the tail itself). We identify the time 
and place of contact with the bend by solving a differential equation similar to (fTl) : assuming the speed of 
the tail is u(t), we know how the bitangent between m and the corresponding ball centered at a changes; in 
particular, we know its length Z(t) as a function of time. The time r* when the head hits the tail is then the 
solution to the equation r* = 1{t*). After solving the equation we know the configuration of the snake at r*, 
and continue moving the mouth to p around the tail. 

The above procedure essentially "develops" the path 7 piece-by-piece. This is consistent with Section [STSl 
where we described how to pull the snake along a given path 7 for the mouth: the pulling was done piece-by- 
piece, which means that it can be performed even if 7 is not given in advance but instead is revealed piece 
after piece. More importantly, using the procedure, we can develop all L-edges incident to mc, piece-by-piece, 
in a BFS manner. We describe this below. 

To initiate the developing, look at the visibility segment m'mc of the L-edge that has led the mouth to 
m,c- The segment is tangent to a slide S 9 m,c- The slide S down from mc (i.e., in the direction consistent 
with m,'mc) is the first (potential) piece of a new L-edge. We extend bitangents from the piece to all other 
1-slides and to all pieces of the spine inflated by 2. These bitangents become the next potential pieces for 
the L-edges. After the bitangents, the next potential pieces are the slides and the spine pieces at which the 
bitangents end. We continue in this way (possibly adjusting the pieces of a particular edge to account for 
snake self-interaction) until for each edge its wiggle segment reaches length L. 

We now bound the time spent on developing all L-edges from mc- Each edge has linear complexity; this 
can be proved identically to Lemma |3.3[ Thus all edges can be grown in polynomial time if the number of 
edges is polynomial. This is what we prove next: 

Lemma 3.5. Let E he the set of all L-edges incident to mc- |IE| = 0{n^). 

Proof. By definition, every edge starts from a path 7 for the mouth of length at most L. Say that paths 
7i , 72 are the same combinatorial type if the sequence of vertices visited by 71 is a subsequence of that for 72 
(or vice versa). Being of the same type is an equivalence relation that splits E into classes. For each class, 
keep only the path with the longest sequence of visited vertices, and identify the class with the path. Let E* 
be the obtained collection of classes. 

Let 7* S E*. Any L-edge that has a path 7 G 7* as its wiggle segment is obtained by extending a 
bitangent from 7 (or equivalently from 7* ) to some slide. Thus the total number of L-edges having a path in 
7* as the wiggle segment is at most the number of bitangents from 7* to the slides, which is 0{n^) since 7* 
is 0(n)-complexity. 

It remains to show that |E*| = 0(1). A standard argument shows that |E*| depends only on the number 
of the holes in P^ reachable by length-L paths from mc (not on the number of vertices): Restrict attention 
to radius-L disk V centered on mc, do vertical trapezoidation of P-'^ n 2?, and let G be the dual graph of the 
trapezoidation. The number of nodes of G of degree higher than 2 is linear in the number of holes (in P^) 
that intersect 2?; call this number H . Transform G to a graph G" by replacing each path in G by an edge 
between degree-3 and higher nodes; the number of nodes in G' is 0{H). Each path from E* is a walk in G' 
that visits every edge of G' at most K — 0(1) times. The number of such walks G' depends only LL , and we 
prove now that LL is constant: 



Claim 1. Let hi . . . h„i be the holes (of P) that intersect D. The Minkowski sum < /ii U . . . U /i„i > has a 
constant number of connected components. 

Proof. If /ii, /i2 are in different connected components of the sum, the distance between the holes is at least 2. 
Within constant distance of mc , there can be only a constant number of points pairwise-separated by distance 
at least 2; thus there is only a constant-size set of holes pairwise in different connected components. D 

This competes the proof of the lemma. D 

3.5 Label! 

We are ready now to traverse the " L- visibility" graph Q of P, searching only the relevant part of Q and not 
building the whole graph explicitly. The label of each node in Q consists of two parts: the distance label 
(storing the distance from s) and the configuration label (storing the snake configuration at which the node 
was reached). That is, a node may have several labels - one per configuration. However, since there is only a 



constant number of different configurations of pulled-taut snakes that may reach the node (Lemma 3.5 and 
Claim fl]), the total number of labels of any node is constant. 

Start from mc = s, and assign distance label to s. The algorithm grows the graph Q whose edges are 
the discovered L-edges and whose nodes are the endpoints of the i-edges. Note that by the definition of 



L-visibility (Definition 3.4), all nodes of G reside on slides. At a generic step, take the node with the smallest 
(temporary) distance label, make the label permanent, and construct L-edges from the node. The endpoints 
of the edges join G and get their (temporary) distance labels and configuration labels. In addition, each 
already existing node over which the mouth passes, gets its distance label updated if the distance label carried 
with the mouth is smaller than the node's current label and the configuration label carried with the mouth is 
the same as the node's configuration label. The search stops when / is reached. 

We now prove that the algorithm terminates in a polynomial number of steps. The visibility segments are 
bitangents between 1-slides and pieces of the wiggle segments. These latter pieces are of two types: (1) slides 
up to layer K and (2) curves that are obtained as the head rolls over the tail, possibly padded by up to 
K layers of the snake. There is in principle an uncountable number of possible pieces of the second type. 
Nevertheless, every step of the algorithm discovers at least one new visibility segment tangent to a piece of 
the first type. Since the number of slides up to layer K is 0(n), there are 0{n^) of such visibility segments. 
Each such segment may be discovered only a constant number of times — once per homotopy type of the 
snake reaching the endpoint of the segment. Thus overall there are O(n^) steps. 

Overall, we have our main positive result: 

Theorem 3.6. Shortest path for a fat hippo can be computed in polynomial time. 

4 Being a Long Snake is Hard 

In this section we prove that if the snake length is not bounded, deciding existence of a path for the 
snake is NP-hard. Specifically, our problem is: Given polygonal domain P and points s, /, find a thick 
non-selfoverlapping s-f path (i.e., a thick s-f wire). 

We show the problem's hardness by a reduction from planar 3SAT. Recall that the graph of a 3SAT 
instance has nodes for all variables and clauses, and two types of edges: (1) a cycle C through all variables, 
and (2) edges connecting every clause to its three variables (Fig. [4]). Planar 3SAT is a restriction of 3SAT to 
instances whose graph is planar; Lichtenstein |18j proved that planar 3SAT is NP-hard. To show the hardness 
of wire routing, we start from an instance / of planar 3SAT; we identify the instance with its (planar) graph, 
and the variables and clauses with the points in the plane into which they are embedded. 




Figure 4: (The graph of) an instance / of 3SAT. The variables are shaded circles, the clauses are hollow 
circles. 




Figure 5: / augmented with parent-child edges, sibling edges, and with variable-clause edges duplicated. 




Figure 6: The closed walks Wout, Wi„; the numbers indicate the order in which edges are traversed by the 
walks. 




Figure 7: The closed walk W; the numbers indicate the order in which edges are traversed by the walk. 
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Augmenting / 

The cycle C splits the plane into two parts; each clause belongs to exactly one part. We say that the clauses, 
edges, etc. inside (resp. outside) C are inner (resp. outer). 

Focus on the outer clauses. Define parent-child relationship between the clauses as follows. The clauses 
that belong to the outer face of the graph I are orphans - they have no parents. Now imagine removing 
an orphan C, together with the edges that connect C to the variables. Any clause C" that now (after C's 
removal) belongs to the outer face of / is a child of C (and C is the parent of C"). Recursively, any clause 
C" is the parent for all clauses that appear on the outer face of / after removal of C" (and all its ancestors). 

We now augment / with new edges. For any parent, the children are angularly sorted around the parent. 
We connect the first and the last child to the parent; if a parent has only one child, connect the parent and 
the child by 2 parallel edges (parallel in the graph-theoretic sense, in the embedding they are not parallel). 
Add also edges between siblings; in particular, connect orphans with a cycle, in the order in which they 
appear on the outer face of / (Fig. [5]) . Finally, add a parallel edge for each clause- variable edge (so that every 
clause is connected to each of its variables with 2 parallel edges) . 

The walks Wout,y^in 

Let Wout be the closed walk that goes through the outer clauses and all variables in the DFS manner, with 
preference to go right (as seen from a clause) and down (i.e., towards C). Specifically, at the "top level", the 
walk contains the cycle through all orphans. In addition, at each orphan C (and in general, at each clause) 
the walk is the DFS traversal of the subtree of C: it goes to the rightmost (as seen from C) variable of C, 
then goes back to C (using the parallel edge), and then - either to the next variable of C or to the rightmost 
child C" of C (whichever is more to the right). At C", the walk recurses down to the rightmost variable of C", 
then goes back to C", and then again - either to the next variable of C" or to the rightmost child of C", etc. 
The walk follows the edge to the sibling or to the parent of a clause C" after all children and variables of C" 
have been visited. In particular, the variables of childless clauses are just visited one-by-one, from right to 
left. Refer to Fig. [e] 

We do an analogous augmentation of the subgraph of / inside the cycle C. Specifically, we choose some 
face F oi I inside C, and let F play the role of the outer face: The orphan clauses are those on the boundary 
of F; the children are defined recursively as the clauses that appear on F after deletion of parents. As before, 
we duplicate variable-clause edges. Let Win be closed walk analogous to Wout- yVin goes through the orphans, 
recursing to variables and children in the DFS manner, with the preference to go left (as seen from a clause) 
and towards C Refer to Fig. [6j 

The walk W 

We now splice the walks Wout, Win, and the cycle C into a single closed walk W through /. Let ab be an 
arbitrary edge of C, and let Cout be an outer clause that belongs to the face of / that has ab on the boundary. 
Remove ab from C, and add edges aCout, bCout (Fig. Uj. Similarly, remove an edge cd from C, and add edges 
cCin, dCin to an inner clause Ci„. 

The walk W starts from following Wout until reaching Cout, at which point it uses the edge Coutb- (This 
may not necessarily be the first time the walk visits Cout', the walk uses Coutb so that the usage is consistent 
with the ordering of edges around Cout - refer to Fig. [7j where Cout is reached first by edge 22, while Coutb is 
only 25th.) From b, the walk follows the cycle C up to c, where it uses the edge cCin to enter inside C From 
Cin, the walk W follows the walk Win all the way around back to Ci„. It then uses Cind to get back to C, 
upon which it traverses C from d to a. At a, the walk uses aCin and follows the rest of Wout from there. 

From 3SAT to wire routing 

As the last step of the reduction, we convert / to an instance of finding a thick wire. For that, we thicken the 
edges of /, turning them into channels of width 2. We replace variables and clauses with gadgets shown in 
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Figure 8: A variable gadget can be traversed in one of the two ways, setting the truth assignment. Depending 
on the way, the path bulges out of the odd (when the variable is set to True) or even (when it is False) 
openings at the top of the gadget, and resp. even/odd openings at the bottom. 
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Figure 9: When a clause gadget is traversed, from left to right, one of the channels leading to variables must 
be used. Otherwise, 3 subpaths go through the top of the gadget leading to a self-overlap. 



Figs. Is] and [9] resp. The connections (channels) between variables and clauses (Fig. [To| ensure that whenever 
a channel from a clause to a variable is used by the wire, the variable satisfies the clause. 

We cut (the channel corresponding to) one of the edges of C, and place the points s, / on the opposite 
sides of the cut. This turns the closed walk W into an s-f path. By our construction, the only s-f wire in the 
instance is one that follows the walk W in /, possibly, omitting some clause- variable edges (channels). Indeed, 
neither channels nor gadgets have any leakage - the wire must follow them through. The only flexibility that 
the wire has is (1) how to traverse the variable gadgets, and (2) which clause-variable channels (not) to use. 
But we know, by clause gadget construction (see Fig.[9|, that at least one channel from every clause must be 
used; moreover, when it is used, the variable must satisfy the clause (see Fig. 10 1. Thus, the s-f wire exists if 
and only if / is satisfiable. 



5 Conclusion 

We showed that for a snake to stay simple, it must grow fat (or else, learn to squeeze). We mention few open 
problems here: 

To minimize snake's squeeze factor, it would be interesting to find a path that minimizes maximum self- 
overlap of a snake. It is easy to see that a snake of thickness 1/2 can follow (without self-overlap) the 
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Figure 10: If a variable does not satisfy a clause the channel between them cannot be used by non- 
selfoverlapping path. 
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shortest (possibly self-overlapping) path for a thickness-1 snake. Is 1/2 best possible? Note that in our 
hardness proof, a non-zero amount of overlap is enforced; what is the largest self-overlap for which the 
problem remains hard? 

Minimizing the path length for one point of the snake is similar to finding "di-optimal" motion of a rod [5]. 
Other objectives are possible: e.g., distance traveled by another point (not the mouth) or the average 
distance traveled (^2-7 ■ ■ ■ 1 doo-optimality, etc.). Does the situation with the snake mimic that for the 
rod: minimizing motion of any point is NP-hard, except for a rod cndpoint? 

What is the hardness of computing paths for a short snake that is not required to be pulled taut? It seems 
that the answer to this question leads to an interesting research direction of "snake packing" : Given a 
polygonal domain, can one layout a length-L snake in it? The problem is NP-hard by a reduction from 
Hamiltonicity of grid graphs; what about, say, simple polygons? 

Can a shortest rectilinear wire be computed efficiently? 
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